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Abstract 

> , 

Preliminary results are presented from the CLEF hydrodynamics simulation, a large 
(N = 2 x 428 3 particles within a 200/i~ 1 Mpc comoving box) simulation of the 
ACDM cosmology that includes both radiative cooling and a simple model for galac- 
tic feedback. Specifically, we focus on the X-ray properties of the simulated clusters 
at z = and demonstrate a reasonable level of agreement between simulated and 
' observed cluster scaling relations. 
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1 Introduction 



As the largest and latest virialised structure s to form, galaxy c lusters are 
especially useful cosmological probes (e.g. see IViana et all 1200.1 and refer- 



ences therein) . Next generation cluster cosmology surveys, such as the XCS 



fjRomer et all 120011 ) . will detect sufficiently large numbers of clusters that 
uncertainties in values of cosmological parameters will be mainly systematic, 
requiring for example an accurate calibration between cluster X-ray temper- 
ature and mass. Such measurements demand an improved understanding of 
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Fig. 1. Left: optical depth image of the gas in a 200 x 200 x 10 h Mpc slice at 
2 = 0. Right: series of zooms showing images of the mass-weighted temperature of 
the gas, from the full box width to an individual cluster. 

cluster physics, therefore realistic numerical simulations of the cluster popu- 
lation are essential. 

In this paper we present a preliminary analysis of the z = cluster population 
within the CLEF hydrodynamics simulation, a large state-of-the-art cosmo- 
logical simulation that, besides gravity and gas dynamics, includes a model 
for the effects of galaxy formation. As we will show, the simulation does a 
reasonably good job at reproducing X-ray scaling relations at z = 0. 



2 The CLEF hydrodynamics simulation 



The CLEF (CLuster Evolution and Formation) hydrodynamics simulation 
(see Fig 1) is a large simulation of structure formation within the ACDM 
cosmology, with the following cosmological parameters: Q m = 0.3, = 0.7, 
VL^h 2 = 0.0238, h = 0.7 and as = 0-9- These val ues are in good agreement 
with recent WMAP analyses ([Spergel et al. 



Initial conditions were generated using a mo dified version of the COSMIC soft 



ware package provided with the hydra code (|Couchman. Thomas and Pearcd. 

1995). The appropriate transfer function, generated using cmbfast ([Seliak and Zaldarriaga . 



1996), was read in and a displacement field generated for a 200 h~ x Mpc co- 
moving box at z — 49. Two regular cubic grids of 428 3 particles, separated by 
half the interparticle distance in each of the x, y and z directions, were then 
perturbed by these displacements to create the initial particle positions. Thus, 
the gas and dark matter particle masses were set to m gas = 1.4 x 10 9 ^ -1 M 
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and mdark = 7.1 x 10 9 h 1 M respectively. 



This initial configuration was then evolved to z — using version 2 of the gad- 
get code (|Springel. Yoshida and Whitel . 120011 ). a hybrid Particle-Mesh/Tree 
gravity solver with a version of Smoothed Particle Hydrodynamics (SPH) 
that explicitly conserves entropy where appropriate. In addition, the gas could 
cool radiatively, assuming a fixed metallicity of Z = 0.3Z Q . Cooled gas, with 
Wh > 10 -3 cm -3 and T < 1.2 x 10 4 K, could either form stars if r > /heat or 
be reheated by stars if r < /heat, where r is a random number drawn for each 
particle from the unit interval and /h ea t = 0.1 is the reheated mass fraction 
parameter. Each reheated gas particle was given a fixed amount of entropy, 
>Sheat = lOOOkeVcm 2 , where 5* = kT/n 2 ^ 3 , which further heats the ICM as 
the particle does work on its surroundings. Further details may be found in 
Kav et all (|2004h . 



3 X-ray scaling relations at z = 



In this paper, we concentrate on comparing a selection of simulated and ob- 
served X-ray cluster scaling relations at z — 0. Clusters were identified by 
first identifying local maxima in the density field and growing spheres around 
these maxima until the average density within each sphere was a fixed factor, 
A, above the critical density, p cr = 3Hq/8tiG. Values of A used will be given 
in each subsection. For the virial density (A ~ 104) there are > 400 clusters 
with kT viT > 1 keV (> 60 above 3 keV). 



3.1 Temperature-mass relation 

We begin by showing in Fig. 2 the relation between hot gas mass-weighted 
temperature (T gas = SjmjTj/Sjmj, where the sum is over all gas particles 
with Ti > 10 5 K) and total mass for a density contrast A = 2500. All clusters 
with M2500 > 3 x 10 14 /i — 1 M are considered. The dashed line is a best-fit 
relation to the clusters for a fixed slope of 2/3, as expected if the clusters form 
a self-similar population. This relation is 

log(fcT gas /keV) = (0.614 ± 0.003) + (2/3) log(M 250 o/M 14 ), (1) 



where M14 = 10 1A h 1 M Q . When the slope is allowed to vary, the best-fit 
relation (solid line) is 

log(A;T gas /keV) = (0.608 ± 0.004) + (0.65 ± 0.01) log(M 250 o/M 14 ), (2) 
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Fig. 2. Gas mass-weighted temperature versus mass, evaluated at A = 2500. The 
dashed line is the best-fit relation with the self-similar slope 2/3. The solid line is 
the best-fit relation, allowing both the normalisation and slope to vary. The solid 
band is the best-fit relation to clusters studied by Allen, Schmidt and Fabian (2001). 

close to the self-similar relation. The thick band is the observed relation de- 



rived by (jAllen. Schmidt and Fabianl . 120011 ). in good agreement with our re- 
sults. 

It is more common in the literature for observed temperature-mass relations 
to be presented at larger radii, using spectroscopic (photon-weighted) temper- 
atures and mass estimates assuming /3-m odel surface brightness profiles and 
polytropic-model temperature profiles (e.g.lNevalainen. Markevitch and Forman 
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2000l: Finoguenov. Reiprich and Bohringeri200lJ : Sanderson et al.ll2003f ). These 
results generally suggest a slope closer to 1/2 than 2/3, attributed to non- 
gravitational processes (see below), and a normalisation that is offset in mass 
from simulation predictions by ~40 per cent. Examining, for example, the X- 
ray emission- weighted temperature- mass relation from our simulation at rsoo, 
Tx — M500 px = T li m i niA(T i )Ti/'E i miniA(Ti), where A is an energy-dependent 
cooling function], we find a similar slope to the observations (0.53) but an off- 
set in normalisation comparable to previous simulations. The cause of this 
offset is likely due to incorrect estimates of cluster masses (e.g. iRasia et al 
2004) and is something we will return to in a future paper. 



3.2 Entropy-temperature relation 

Galaxy formation increases the entropy of intracluster gas, producing a rela- 
tionship with temperature that is flatter than the self-similar scaling (S ocT). 
We plot this relation in Fig. 3, again using an X-ray emission- weighted temper- 
ature for each cluster. Two radii are considered (0.1r 2 oo an d ^500) an d only clus- 
ters with kTx > IkeV are studied. Again, t he simulated clusters are in reason- 
ably good agreement with the observations (jPonman. Sanderson and Finoguenov . 
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Fig. 3. Entropy versus X-ray emission-weighted temperature at 0.1r2oo (lower 
points) and rsoo (upper points), with the solid lines being fits to these data. Crosses 
are data from Ponman, Sanderson and Finoguenov (2003) and dashed lines are 
self-similar scalings, normalised to their hottest clusters. 

2003j ). containing an excess of entropy that is larger in smaller systems. For 



0.1r 200 

log(Sy keV cm 2 ) = (2.14 ± 0.008) + (0.46 ± 0.03) log(A;T x / keV) (3) 
and for r 500 

log(S/keVcm 2 ) = (2.84 ± 0.003) + (0.63 ± 0.01) log(A;T x / keV). (4) 



3.3 Luminosity-temperature relation 

Finally, we show bolometric X-ray luminosity versus X-ray emission-weighted 
temperature in Fig. 4. Again, only clusters with kTx > IkeV are c onsid- 
ere d. Symbols with errors are observational data from Markcvit 3 f)l998h 



and Arnaud and Evrardl (1999). To remain approximately consistent with this 



data, emission from within 50 h 1 kpc of our simulated cluster centres is omit- 
ted. 

The dashed line is a best-fit relation for a fixed slope equal to 2 (self-similar), 
clearly a poor fit to the observations. When the slope is allowed to vary, the 
best-fit relation (solid line) is 

log(L x /L 40 ) = (1.89 ±0.01) + (3.84 ± 0.05) log(A;T x /keV), (5) 



where L 40 = 10 40 h 2 ergs 1 , considerably steeper than the self-similar case. 
In fact, the simulated relation is not adequately described by a power law since 
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Fig. 4. Bolometric X-ray luminosity versus emission-weighted temperature in the 
1.5-8.0 keV band. Symbols with error bars are data from Markevitch (1998) and 
Arnaud and Evrard (1999). The dashed line is the best-fit relation for a fixed slope 
of 2. The solid line is the best-fit relation for kTx > IkeV and the dot-dashed line 
for kT x > 3keV. 

the local gradient becomes progressively natter with increasing temperature. 
Fitting clusters with kTx > 3keV (dot-dashed line) yields 

log(L x /L 40 ) = (2.47 ± 0.03) + (2.88 ± 0.05) log(A;T x / keV), (6) 



in reasonable agreement with the observations a lthough the norm alisation is 
a bit too high. Better agreement was reached bv lKav et al.l (j2004h . who used 
a slightly smaller gas fraction (0.15 rather than 0.162 used here, which is a 
closer match to the WMAP value). It is likely however that fine tuning of 
the feedback model parameters would improve the agreement between the 
simulated and observed relations. 
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